Thouless energy of a superconductor from non local conductance fluctuations 
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We show that a spin-up electron from a normal metal entering a superconductor propagates as 
a composite object consisting of a spin-down hole and a pair in the condensate. This leads to a 
factorization of the non local conductance as two local Andreev reflections at both interfaces and 
one propagation in the superconductor, which is tested numerically within a one dimensional toy 
model of reflectionless tunneling. Small area junctions are characterized by non local conductance 
fluctuations. A treatment ignoring weak localization leads to a Thouless energy inverse proportional 
to the sample size, as observed in the numerical simulations. We show that weak localization can 
have a strong effect, and leads to a coupling between evanescent quasiparticles and the condensate 
by Andreev reflections "internal" to the superconductor. 

PACS numbers: 74.50.-|-r,74.78.Na,74.78.Fk 



I. INTRODUCTION 

Correlated pairs of electrons can be manipulated in 
solid state devices by extracting Andreev pairs from 
a conventional superconductor, being a condensate of 
Cooper pairs. This process is known as Andreev 
reflection^ at a normal metal / superconductor (NS) in- 
terface. One considers the future realization of devices 
designed for manipulating separately one of the two elec- 
trons of an Andreev pair and see the feedback on the 
other electronSiSi^iSiSii. The question arises of exploring 
experimentally and understanding theoretically the prop- 
erties of the simplest of these devices: a source of spa- 
tially separated Andreev pairs propagating in different 
electrodes forming, in short, "non local" Andreev pairs. 
The possibility of realizing a source of non local Andreev 

pairs has indeed aroused a considerable interest recently, 
both theoretical'r'i4r^"i7iVd"diii2,i.3,i4,i.5,i6,i7,i8,i9 ^^^ 

experimentalSSiSi. 

In a theoretical prediction prior to the 
experimentsSSiSi, Falci et alm^ have obtained from 
lowest order perturbation theory in the tunnel ampli- 
tudes a vanishingly small non local signal with normal 
metals. Russo et al?^ have obtained on the contrary a 
sizeable experimental non local signal in a three terminal 
structure consisting of a normal metal / insulator / 
superconductor / insulator / normal metal (NISIN) tri- 
layer. The goal of our article is to provide a theory that, 
together with Ref. [I^ contributes to the understanding 
of this experimenliSi, as well as related possible future 
experiments on non local conductance fluctuations, and 
be consistent with the other available experiment by 
Beckmann et al. with ferromagnetsSS. 

Falci et ali^ have discussed the two competing chan- 
nels contributing to non local transport. An incoming 
electron in electrode "b" can be transmitted as an elec- 
tron in electrode "a" , corresponding to normal transmis- 
sion in the electron-electron channel (see the device on 
Fig. n for the labels "a" and "b" ) . Conversely, it can be 
transmitted as a hole in electrode "a" while a Cooper pair 
is transfered in the superconductor. The later transmis- 



sion in the electron-hole channel corresponds to a dom- 
inant "non local" Andreev reflection channel that can 
lead to spatially separated, spin entangled pairs of elec- 
trons. The outgoing particles in the two transmission 
channels have an opposite charge, resulting in a differ- 
ent sign in the contribution to the current in electrode 
"a" . With normal metals, not only have the elastic co- 
tunneling and crossed Andreev reflection an opposite sign 
in the non local conductance, but they are exactly oppo- 
site within lowest order perturbation theory in the tunnel 
amplitudes. 

It was already established that non local trans- 
port is dominated by elastic cotunneling for localized 
interfacesi&. The superconductor can essentially be re- 
placed by an insulator for a very thin superconductor 
connected by tunnel contacts to a normal metal (assum- 
ing that the superconductor can still be described by 
BCS theory). We show that this picture breaks down 
if the superconductor thickness is larger than the coher- 
ence length because transport is mediated by composite 
objects made of evanescent quasiparticles and pairs in 
the condensate. 

On the other hand, we find that small area junctions 
are controlled by a different physics with fluctuations 
of the non local conductance. We find on the basis of 
an evaluation of the diffuson in a superconductor, that 
the Thouless energy is inverse proportional to the system 
size, which matches our numerical simulations. We find 
also a possible large coupling to the condensate provided 
by weak localization in the superconductor. 

The article is organized as follows. The factorization of 
non local processes as two local Andreev reflections and 
a non local propagation is discussed in SecEJ The fac- 
torization of the non local conductance is illustrated in 
Sec, mil in the case of one dimensional models (the Blon- 
der, Thinkham, Klapwijk (BTK) model^ and a Green's 
function model) . The Thouless energy of non local con- 
ductance fluctuations is examined in Sec. II VI on the basis 
of the evaluation of the diffuson. Numerical simulations 
are presented in Sec. The role of weak localization is 
pointed out in Sec. IVII Concluding remarks are given in 
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FIG. 1: (Color online.) Schematic representation of the elec- 
trical circuit corresponding to the NISIN double interface in- 
terpolating between a localized contact for L ~ Xf (with 
Xf the Fermi wave-length) and and extended interface for 
L S> Lth{u!), where Lth is the Thouless length correspond- 
ing to the energy eVt. The current la through electrode 
"a" is determined in response to a voltage Vt on electrode 
"b", with Va = 0. The dimensions R (superconductor thick- 
ness) and L (dimension of the junction) are shown on the 
figure. The available experiment by Russo et aim— corre- 
sponds to extended interface with L of the order of O.S/xm 
and R~15^ 200nm. 
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II. FACTORIZATION OF THE NON LOCAL 
RESISTANCE 

A. Existing results for eVJ, > Eth{L) 

The diagram corresponding to the non vanishing low- 
est order process of order T"* (with T the normal trans- 
parency) is shown on Fig.Oi. This diagram is local with 
respect to excursions parallel to the interfaces if the bias 
voltage energy eVb is much larger than the Thouless en- 
ergy Eth{L) associated to the dimension L of the junction 
parallel to the interface (see Fig. QJ, as it is the case in 
the experiment by Russo et alm^. The corresponding non 
local conductance 



— dl 



(1) 



where la and Vb are defined on Fig. Q] is given bjsiSiiSi 



SaAVb)=-jN,,.T ^^^^^^_^^^^^^ 



exp 



d 



(2) 

where Nch is the number of conduction channels, A the 
superconducting gap, ^ the superconducting coherence 

(S) 

length, le the superconductor elastic mean free path. 



T the normal local transparency, and the overline is an 
average over disorder. The local Andreev conductance is 
given by 

SUV.)^2'^N.,T^^^-^, (3) 

where we used the ballistic result without disorder be- 
cause of the condition eVb » Eth [L) . 

The resistance matrix probed in the experiment^i is 
the inverse of the conductance matrix calculated theo- 
retically: 



■Ra.aiVb) na,b{Vb) 

7^b,a(l4) 7^b,6(Vfc) 



Q_a,a{yb) Q_a.b{Vb) 

GbAVb) QbAVb) 



-I -1 



(4) 



from what we deduce that the non local resistance 
7?.Q,b(Vfc) is given by 



na,b{Vb) 



-Qa,b{Vb, 



[QxociVbW -Qa.b{Vb)gbAyb)' 



that simplifies into 

Tla.biVb) 



GaAVb) 
[GlociVbt 



(5) 



(6) 



if the thickness of the superconductor is larger than the 
superconducting coherence length ^, and leads to 



na,b{Vb, 



1 h 

A^ - {eVbf 
A2 



exp(-2i?/e). (7) 



The non local resistance at low bias is positi ve ( domi- 
nated by elastic cotunneling), as found in Ref. Ila The 
case of extended interfaces is addressed in Ref. 119. 



B. Case eVb < Ea,{L) 

Now, if the bias voltage energy eVb is smaller than 
the Thouless energy Eth{L), the lowest order diagram of 
order T'^ becomes non local in the normal electrodes (see 
Fig. 135). The diagram on Fig. I^b corresponds to two 
Andreev reflections in the normal electrodes, connected 
by a propagation in the superconductor, so that the non 
local conductance factorizes into 



Qa,b{Vb) 



S{Vb) [Qioc{Vb)f 



N, 



ch 



(8) 



where S{Vb) is a transmission coefficient of the supercon- 
ductor. Using Eq. (jHJl, we find that the crossed resistance 
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FIG. 2: (Color online.) The diagrams representing the low- 
est order processes of order T**. The diagram (a), local with 
respect to a prop agation in the normal electrodes, was intro- 
duced in Ref.lld The diagram (b) is its non local counterpart. 
"1" and "2" correspond to the electron and hole Nambu la- 
bels. The diagram on (b) factorizes in two Andreev reflections 
at both interfaces, and a non local propagation in the super- 
conductor. The electron line crosses 8 times the interfaces, so 
that the diagrams are of order T*, where T oc {t/eFY is the 
normal transparency, with i the tunnel amplitude and ep the 
Fermi energy. 



does not depend on the local conductances. The scaling 
between the local and non local conductances is tested in 
Sec, mil for the generalization of the model of reflection- 
less tunneling at a single interface introduced by Melsen 
and Beenakker— . 



The factorization of the Andreev reflections at both 
interfaces suggests that part of the current is carried by 
pairs in the condensate. We thus arrive at the notion of 
the transport of a composite object made of an evanes- 
cent quasiparticle and a pair in the condensate: an elec- 
tron from a normal electrode is transmitted in the super- 
conductor as a quasi-hole and a pair in the condensate 
(see Fig. E^). The consequences of this qualitative pic- 
ture are considered below. 



Finally, we note that the factorization of two Andreev 
reflections at the interfaces is also valid if eVf, ^ Et\^{L) 
(see Sec. Ill A|l . as in the experiment by Russo et alm^. 
This is because the normal Green's functions are vanish- 
ingly small at zero energy in a superconductor. 



III. ONE DIMENSIONAL MODELS 
A. Blonder, Tinkham, Klapwijk (BTK) approach 

1. Non local conductance 

Let us consider now a one dimensional model of 
NISIN double interface within the BTK approach^"* (see 
Fig. Eli). The goal is two- fold: i) obtain the expres- 
sion of the pair current in the superconductor, and ii) 
test the factorization of the non local conductance in the 
case of the model of reflectionless tunneling introduced 
by Melsen and Beenakker--^. 

The gap of the superconductor is supposed to have a 
step-function variation as a function of the coordinate z 
along the chain: A(z) = A9{z + R/2)0{R/2 - z), and 
we suppose (5-function scattering potentials at the inter- 
faces: V{z) = HS{z + R/2) + HS{z - R/2]M.. The inter- 
face transparencies are characterized by the parameter 
Z = 2mH/h^kp, where vp = hkp/m is the Fermi ve- 
locity, with m the electron mass and kp the Fermi wave- 
vector. The one dimensional model is a simplified version 
of the genuine three terminal geometry with a supercur- 
rent flow. The current in the normal electrode "a" is not 
equal to the injected current in electrode "b" because 
part of the injected current has been converted in a su- 
percurrent. 

The unknown coefficients in the expression of the wave- 
function are determined from the matching conditions at 
the interfaces?? . Of particular interest are the amplitudes 
a'{kpR) and b'{kpR) of transmission in the electron-hole 
and electron-electron channels from one normal metal to 
the other, corresponding respectively to elastic cotunnel- 
ing and non local Andreev refiection. Assuming R^ ^, 
we expand a' and b' to lowest order in exp(— _R/^), to 
find the transmission coefficients 



-— \a{kpR)\ = 

ZTT 

[ 2Z* 2Z6 + 2Z» ^ 



^"^ d{kpR).,,., ^,,2 



(10) 



-2i^/? + o (e-4^/«) 



■\b'{kpR)\^ 



(11) 



27r 

1 1 5 

2Z* ^ 2Z^ ^ 4Z8 

at w = 0. We deduce the first non vanishing term in the 
large- i?, large-Z expansion of the non local transmission: 






T' 



d{kpR) 



4Z8 



27r 



{\a'{kpR)\'-\h'{kpR)\') (12) 



O 



(e-4fl/c) . 



In agreement with the Green's function approacliiSii^ 
corresponding to the diagrams on Fig. [3 the non local 
conductance is dominated by elastic cotunneling and ap- 
pears at order Z~^ ^ T^. In agreement with Ref. Ha we 
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FIG. 3: (Color online.) Schematic representation of (a): the formation of composite objects at the interfaces with a ballistic 
superconductor; (b): the sequential conversion of the composite objects in the bulk of the superconductor in the absence of 
weak localization; (c): the scattering induced by weak localization. 
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FIG. 4: (Color online.) Schematic representations of the one 
dimensional models: (a) the BTK model for NISIN and (b) 
NINISININ junctions, and (c) the tight-binding model for 
NISIN and (d) NINISININ junctions. 



find no non local Andreev reflection for highly transpar- 
ent interfaces corresponding to ^ = 0. 



B. Reflectionless tunneling 

1. BTK approach 

To discuss the form @ of the crossed resistance, we 
include now multiple scattering in the normal electrodes 
and consider two additional scatterers at positions zi ~ 
—Li/2 in the left electrode and Z2 — ^2/2 in the right 
electrode, described by the potentials V'{z) = H'S{z — 
zi) + H'5{z — Z2), and leading to the barrier parameter 
Z' = 2mH' /h^kp (see Fig.0J) for the definitions of Z and 
Z'). This constitutes, for a double interface, the analog 
of the model introduced by Melsen and Beenakker— for 
a single interface. We average numerically the non local 
transmission coefficient over the Fermi oscillation phases 
(y3i — kp{R — ii)/2, if = kpR and (^2 = kp{L2 — R). 

The variations of the non local conductance at zero 
bias as a function of Z' for a fixed Z are shown on Fig. [SJ 
as well as the corresponding non local conductance for the 
NINISIN junction. The negative non local conductance 
at small Zi for the NINISININ junction disappears when 
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FIG. 5: (Color online.) Variation of the non local conductance Ga.b (in units of e^ /h) for the junction on Fig. |1J:, with 
Z' = Zi — Z4, and Z = Z2 — Z^ — 10. (a) ... (f) correspond to an increasing values of the precision in the evaluation of 
the Fermi phase factors related to the superconductor. We have also shown the much smaller non local conductance of the 
NINISIN junction, as a function of Zi for the NIN contact, with the same value of Z for the NIS contacts. 
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Green's functions: scaling between the local and non 
local conductances 
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FIG. 6: (Color online.) Correlation between the nonlocal 
conductance (a;-axis) and the local conductance (j/-ajds) for 
the peaks similar to Fig. I^as a function of t' /T for a fixed 
t/T = 0.05 and (a) R/^ = 5, (b) R/^ = 4, (c) R/^ = 3, 
and (d) 7^/^ = 2. The solid line is a fit to Qa,b{t /T) ex 
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increasing the precision of the integrals. The variation of 
the non local conductance on Fig. shows a strong en- 
hancement by the additional scatterers, like reflect ionless 
tunneling at a single NIS interfaceSS,. 



Considering the tight-binding model within Green's 
functions, the variation of the non local conductance of 
the NINISININ junction as a function of t' for a fixed 
t (see Fig. 01i) is similar to the BTK model. Imposing 
the same normal conductance in the BTK and in the 
tight-binding models leads to Z = (1 - {t/Tf)/{2t/T), 
where T is the bulk hopping amplitude. The identifi- 
cation of Z to t/T results in a good (but not perfect) 
agreement for the non local conductance when the tight- 
binding and BTK results are rescaled on each other. The 
non local conductance Ga.b{Vb = 0, t/T, t' /T) is shown on 
Fig. as a function of the local conductance Q\oc{Vb = 
0,t/T,t'/T), fitted by ga,b{Vb = 0,t/T,t'/T) = S{Vb = 
0) [^ioc(Vh = 0,t/T,t'/T)]^, corresponding to Eq. ® for 
A'ch = 1- The scaling is very well obeyed, showing the 
validity of form Q of the crossed resistance involving the 
destruction of a pair in the condensate at one interface, 
its propagation in the superconductor and its creation at 
the other interface. 











Oa 






















■ 








^ 


M 


^ 


- 






























' 


' 










' 








- ^pJ 


«R 









































FIG. 7: (Color online.) Schematic representation of the recur- 
sive calculation of the matrix diffuson in the ladder approxi- 
mation used by Smith and AmbegaokajS^. The integrals over 
the wave- vectors are carried out recursively, starting from the 
right of the diagram. The two horizontal black lines corre- 
spond to the advanced and retarded Green's functions. The 
green dots correspond to the disorder vertices, and the verti- 
cal red lines correspond to the impurity lines. We have shown 
the Nambu labels {ua, o-r) and {(j'a, cth) at the extremities of 
the diffuson. 



IV. THOULESS ENERGY OF A DISORDERED 
SUPERCONDUCTOR 

A. Relevance to experiments 

We consider now non local conductance fluctuations. 
The total non local transmission coefflcient is given by 
Ttot — Te-e — Te-h, where T^-e and T^-h are the trans- 
mission coeSicients in the electron-electron and electron- 
hole channels respectively. As discussed in Sec. m one 
has Ttot = but 



(Ttot) = {Te- 



.f + {Te-hf 



2{Te-eTe-h) 



(13) 



one less rung, but with a different 2x2 matrix at the ex- 
tremity. The four parameter recursion relations reduce 
to a matrix geometric series in the sector (fo,Ti), and 
to another matrix geometric series in the sector {t2,t^)^ 
where f„ are the four Pauli matrices, with tq the identity. 
More precisely, we define the four matrix diffusons 



^ q^Su V n 



^2^ 



f3G(k, a;)f„G(k -H q, cj + (5a))f3, 



(14) 

where n = Q, ...,3, q — |q| is the modulus of the wave- 
vector, and Suj is small compared to the energy lo. The 
microscopic disorder scattering potential is given by v^ = 
ATTtp/Te, with ep the Fermi energy and Te the elastic 
scattering time, related to the elastic scattering length Zg 
by the relation l^ ~ vpTe- We find 



T^q.SujiTl) = X (ujAtq - UJ^fl) 

in the sector (fg , ti ) , and 



(15) 
(16) 



V 



q,6L 



,{~tT2)=X 



SVqASujVA^ 
2^ 



-T3- 



(17) 



^iA' ~lo')t,) 

'i>q,Su>{T''i) ^ XVqSuJ 



VlA6uj 
44 



-T-3 



^^ 



Up-T-l 



in the sector (j2, T3). We used the notation 



1 



3I?o(A2 



VA^ 



V^q^ 



loSlu 



2x/A2 



where Vq is the diffusion constant. 



- UJ 

(18) 



Inspecting the corresponding 



lowest order diagrams 

shows that T^-eTe-h = —{Te-e) , where we suppose that 
the normal metal phase coherence length is vanishingly 
small, therefore avoiding the specific effects of extended 
interfacesiS. The root mean square of the non local con- 
ductance fluctuations is thus proportional to {e'^/h)T^ 
while the average non local conductance is proportional 
to {e-^ / h)T'^ Nch (see Eq. EJ. The fluctuations are impor- 
tant for small junctions such that T'^Nch ^ 1- 



B. Diffusons in a superconductor 

1. Evaluation of the diffusons 

Let US first evaluate the Thouless energy of a super- 
conductor in the absence of crossings between diffusons. 
Smith and Ambegaokar— ' start from one extremity of 
the ladder diagram and calculate recursively the integrals 
over the wave-vectors. Once the right-most integral on 
Fig.[7|has been evaluated, one is left with a "ladder" with 



2. Non local transmission coefficient 

The relation between the diffusons in the superconduc- 
tor and non local transport is provided by the non local 
conductance ^. The non local conductance ^^^(i^) of 
order T^ is related to the transmission coefficients ac- 
cording to 

^2 



t(2)/ 



(1,1), 



-^(2,2) 



Ga,U-)-^ ^a.i/H-ni.i/M 



with 



T, 



(<^'a^'^'r')(,.,\ _ 



2,2 
F 



, ■■'iuo) = T'e 

(aA.iyn) ^ I 



d\ 

(27r)3 



,iq.R 



(19) 



(20) 



(27r)= 



C^^'^'a (k, u;)G""«'^«(k + q, uj). 



The notation T,^''^'^-!(uj) corresponds to the transmis- 
sion coefficient related to a diffuson with the Nambu la- 
bels {<7A,c7ii) for the advanced and retarded propaga- 
tors at one extremity, and the Nambu labels (cr^, ct^) at 



the other extremity (see Fig.[7Il. The transmission coef- is evaluated by discarding the Nambu components of 



(1 1) (2 2) 

ficients T,i\\{iLj) and T,i\\{lj) encode elastic cotunnel- 
ing and non local Andreev reflection respectively. With 



the notations in Sec. IIV A1 we have Tri\\{oj) = T^-ei^) 
for transmission in the electron-electron channel, and 

(2.2) 

T/i \-.{uj) = Te-hii^) for transmission in the electron-hole 

channel. We deduce from Eq. l(TH|) that Ga bi'^) = 0- the 
average non local conductance vanishes to order T^, in 
agreement both with Sec. E] and with an early worki^ 
in the disordered case. The transmission coefficients 

(1.2) (2,1) 

T/i 2)i'^) ^'^d T,2i){uj) involve the propagation of a pair 
in the condensate in parallel to the quasiparticle chan- 
nels, as in the diagram on Fig. [2l3. 



C. Thouless energy 

The Thouless energy is defined from the non local con- 
ductance fluctuations by the decay of the autocorrelation 
of the non local conductance 



( Ga,b{(^)Ga.bi^ + ^^) - Ga.bi^)Sa,bi^ + ^^) )" (21) 



as Slo increases, where {...)uj denotes an average over the 
energy a; in a given window. 

The autocorrelation of the non local conductance de- 
fined by Eq. (|21|) is related to the autocorrelation of the 
transmission coefficients 



^ {a-i,a-2) ^ ' (0-3, 0-4) V / 

-^(^l,-T2)(^)^(rT3,a4) (^ + ^'^ri- 



(22) 



More precisely, the non local conductance to lowest 
order in the tunnel amplitudes is given by 

n I \ /I r 1-1.^ \.\,R 1.2, A 2,l,fll /f,o\ 

GaA^)^^\9a.b 9b,a 'da.b 9b'a J' (23) 

where "1" and "2" refer to the electron and hole Nambu 
components, "A" and "R" stand for advanced and re- 
tarded, g^\ is a propagation from "a" to "b" in the 

electron-electron channel, and g^'^ in the electron-hole 
channel. The prefactor A, not directly relevant to our 
discussion, can be found in Ref. [l^ We find easily 



Ga.biuj)Ga.b{^ + S^) = -^ 



_ /|2 Y^ i{ki-k2+k3-k4)B. 



ki ,...,^4 



(24) 





gtii^^i ^ ^) 9i, iik2, uj) 9i. liks, uj + 6u})g^^^{k4,uj + Suj) 


— ( 

— ( 


]^^^{ki,u;)g^^i{k2,uj)gt^2ik3,^ + Suj)g^^^{k4,uj + 6lj) 


j^2{ki,uj)g2,i{k2,uj)gt,iih,i^ + Suj)g^j^{k4, uj + Suj) 



+ 9i.2iki,^)g2,iik2,uj)gf2ik3,^ + 5uj)g^'-^{k4, u + 5uj) 
The quantity 



Sa,b{'^,LO + 5uj) = Ga.bi^)Ga.bi^ + S^) - Ga,b{^)Ga,b{^) 

(25) 



the type g^^{ki,LL!)g^2{k2,'-^ + Su), much smaller than 

5M(fcii ^)5m(^2, w -t- Su) and gf2{ki,uj)g2,i{k2,(^ + Suj) 
if Lu is small compared to A. In addition, we use 



Ga,bi^)SaA^) — SaA^ + SLj)Ga,bi^ + ^w) (26) 

within the small energy window that we consider. We 
obtain 



Sa,b{^,UJ + Suj) 



(27) 



(r(^^5 {R,u;,Su;)y + 2 (r^ff ^ (i?, ^, ^^)) ' 
(tII:1^^ {R,u;,0)y~2 (rg'5 (i?, a;, 0)) ' , 



where the Fourier transforms with respect to the spatial 



variable of T^ 'J {R, uj, Suj) and T^^ 'J (i?, uj,Suj) are given 



by 



^(M)' (9' ^,Suj) ^ T'^^f^ {q, uj, 5uj) 



(28) 



3P2(A2_^2)[^A2 



Voq^/A - ujSuj/2VA^ - cj2] 



deduced from Scc. IIVBTI The notation R stands for the 
distance between the contacts "a" and "b". Taking the 
Fourier transform of Eq. (|28|) , we obtain 



(29) 



Sa,b{uJ,Suj) 


Ga,b{^)QaA^) 


( f-R ujSuj \ 
= exp iv3— -— -p5 ^ - 


dephases above the Thouless energy 


Ec - Sujc - 


St: A^-uj^ ^ 

y/3 UJ R' 



(30) 



for energies uj large compared to Ay^S^/R, so that Ec is 
much smaller than uj. 



V. NUMERICAL RESULTS 

A. The different length scales in the simulations 

The non local transport simulations are carried out in 
a quasi-lD geometry, on a strip of longitudinal dimen- 
sion L and of transverse dimension Ly, corresponding to 
M transverse modes. We calculate non local transport 
along the z direction. The trilayer geometry with an as- 
pect ratio similar to the experiment by Russo et alm^, 
would require much larger system sizes to have a reason- 
able separation between the different length scales in the 
y direction while L is much larger than Ly. The relevant 
length scales in the diffusive regime are given by increas- 
ing order by the Fermi wave-length \p, the elastic mean 
free path l^, the superconducting coherence length ^ and 
the sample size. 
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FIG. 8; (Color online.) Energy dependence of the supercon- 
ducting transmission coefficient T'(a;/A) [defined by Eq. 1)31^ ] 
through a diffusive superconductor on a strip with Ly/ao = 10 
and L/ao = 100, in the limit of small disorder. The bold line 
corresponds to the ballistic result, and the two other traces 
correspond to two realizations of disorder with W/to — 0.5 
and le/ao — 500. The ballistic coherence length is ^/ao — 33. 
This figure has been obtained with a method based on the 
inversion of the Dyson matrix. 
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FIG. 9: (Color online.) Energy dependence of the supercon- 
ducting transmission coefficient T'{u!/A) [defined by Eq. 1)31^ ] 
through a diffusive superconductor on a strip with Ly/ao = 10 
and L/ao = 100, in the diffusive limit. The two traces cor- 
respond to different realizations of disorder corresponding to 
W/to = 3 and l^/ao — 16. This figure has been obtained with 
a method based on the inversion of the Dyson matrix. 



C. Thouless energy of a disordered superconductor 



B. Ballistic system and small disorder 

We use typically M = 10, L/qq ranging from 80 to 100 
for a methods^ based on the inversion of the Dyson ma- 
trix, and much higher values of L/gq (in units of the tight 
binding model lattice spacing ao), for a complementary 
method consisting in connecting together several conduc- 
tors by a hopping self-energy, given the Green's functions 
of each conductor evaluated by the inversion of the Dyson 
matrix. Disorder is introduced as in the Anderson model 
by a random potential between —W and W. 

The normalized transmission coefficient T'{uj) that we 
calculate numerically is related to the non local conduc- 
tance by the relation 



r'(^) 



h 



T-'s!^li^) 



(31) 



where G„ l{uj) 



^ a^h ^ 



is the contribution of order T^ to the 
non local conductance, with T the normal transparency. 
The transmission coefficient T'{uj) defined by Eq. H31(l 
fluctuates around zero as a function of energy because 
the wave- vectors of the different channels vary with en- 
ergy. The characteristic energy scale in the oscillations of 
the transmission coefficient is the ballistic normal state 
Thouless energy associated to the dimension L (see Fig.|Sl 
in the forthcoming section |VC|I . 



Figs. IHl and |51 show the energy dependence of the su- 
perconducting transmission coefficient T'(lu/A) defined 
by Eq. H31|l . The fiuctuations of the transmission coef- 
ficient are close to the ballistic limit result in the limit 
of small disorder (see Fig. IS)). We obtain regular fluctua- 
tions of the transmission coefficient of a superconductor 
in the diffusive limit where the normal transmission co- 
efficient is characterized by fluctuations (see Fig.EJ). We 
used a large number of realizations of disorder at a single 
energy to show that the transmission coefficient averages 
to zero because of disorder. This shows that the regular 
fluctuations in the disordered system are genuinely re- 
lated to disorder, and do not have the same origin as in 
the ballistic system. 

To characterize the regular fluctuations, we calculate 
the normalized autocorrelation of the transmission coef- 
ficient C{6iLj) — {A{uj,5lo)/B{u),5u)))^, with 



A{lo,5lo) = T'{uo + 5uj)T'{lo) 
B{uj,Suj) ^ Jt'{uj + Suj)T'{uj), 



(32) 
(33) 



where {...)i^ is an average over the energy w, as in 
Eq. H21|l. The autocorrelation C{Suj) is characterized by 
oscillatory damped oscillations (see Fig. I10|) , in contrast 
to the autocorrelation of conductance fluctuations in the 
normal case that is damped without oscillations. The en- 
ergy scales Ec,i and £^£,,2 related to period of oscillations 
and to the damping increase as the system size decreases 
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FIG. 10: (Color online.) Autocorrelation function of the 
transmission coefficient [see Eqs. H32|l -I)33| | ]. for {W/to — 3, 
L/ao = 80), [W/to = 3, L/ao = 100), and {W/to = 4,1/ ao = 
100). W/to — 3 corresponds to le/ao ~ 16, and W/to — 4 
corresponds to h/ao — 10. The errorbars are smaller than 
the size of the symbols. 
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FIG. 11: (Color online.) Variation of the Thouless energy 
Ec,i/A as a function of the system size L/ao in a log-log 
plot, for (a): W/t = 1 and L/ao ~ 127 (b): W/t = 2 and 
le/ao — 34. The solid lines are a fit to Ec,a/A ^ 1/L. 



(see Fig. I10() . in agreement with the expected behavior 
for Thouless-Hke energies. Going to larger system sizes, 
we find that i?c,i scales like the inverse of the sample size 
(seeFig.inj- 

The comparison between Fig.|51for {W/to — 3, le/ao ~ 
16), and similar data for {W/to = 2, le/ao — 34) and 



{W/to = 4, le/ao — 10) show that J5c,i and Ec,2 have 
a weaker dependence on the elastic mean free path than 
for a normal diffusive system. 



VI. EFFECTIVE SCATTERING INDUCED BY 
WEAK LOCALIZATION 

We find a formal analogy between the calculation of 
the non local conductance fluctuations in the preced- 
ing section, and the linear response theory of collective 
j^Q^gg27,28j29^ Namely, the non local conductance fluc- 
tuations can be viewed as a generalized susceptibility in 
linear response, but otherwise the two models involve 
rather different physical quantities. We show now that 
weak localization can induce additional couplings to the 
condensate. Evaluating all the Nambu labels at the two 
three-diffuson vertices (see Fig. IT^ ) is possible in the 
limit w <C A because of the constraint that the normal 
local Green's functions can be discarded in this limit that 
corresponds "internal" Andreev reflection processes as on 
Fig. 121:. The diagrams on Fig. 112b then defines a set of 
transmission coefficients modified by weak localization. 



With the notations A = T', 



(1,2) 



JJ _ y(l.l) 
~ ^(1,1) ' 

coefficients 



A 

B 

C 
D 



(1,2) ■ 



B = T, 



(2,1) 



(1,2) 



c = r; 



(2,2) 



(1,1)^ 



we find for the "renormalized" transmission 

A + X[AB{A + Bf + D^{A^+B^) (34) 

2ABC^] + 0{X'^) 

B + X[AB{A + B'f + C^{A^ + B^) (35) 

2ABD^] + ©(A^) 

C + XD{A + B){D^+3C^) + 0{X'^) (36) 

D + XC{A + B){C^+3D^) + 0{X^), (37) 



with the perturbative parameter A ~ (r^ei^A)^* that can 
turn out to be large. Weak localization can thus lead to a 
large effective scattering for the processes on Fig. ^ with 
multiple imbricate Andreev reflections providing a cou- 
pling between the condensate and the evanescent quasi- 
particle channels. 



VII. CONCLUSIONS 

To conclude, we have provided a theory of non local 
conductance fluctuations at normal metal / superconduc- 
tor double interfaces. First, reconsidering the case of the 
average non local conductance, we confirm that the cen- 
tral role is played by higher order processes in the tunnel 
amplitude. We found that for these processes part of 
the non local current circulates as pairs in the conden- 
sate, not only as evanescent quasiparticles. The crossed 
conductance at zero bias factorizes in the Andreev con- 
ductances at the two interfaces, and a factor related to 
the propagation in the superconductor. This factoriza- 
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FIG. 12: (Color online.) (a) One of the Gorkov-Larkin-Khmelnitskii-Hikami (GLKH) boxes 



I in the superconducting case, 



dressed by a diffuson. (b) The self-crossing of a diffuson with the superconducting GLKH box. (c) Another representation of 
(b), with two additional diffusons at the entry and exit of the GLKH box. "A" and "R" stand for advanced and retarded. 



tion was tested in the context of a model of reflectionless 
tunneling^'^ . 

On the other hand we found numericaUy regular fluc- 
tuations of the non local conductance. The Thouless en- 
ergy inverse proportional to the system size obtained in 
the simulations can be interpreted by a model ignoring 
weak localization. Alternatively, an energy scale inverse 
proportional to the system size could have received an 
interpretation in terms of the Anderson3?-Bogoliubo\i2£ 
collective mode, which is not in contradiction with the 
fact that weak localization can induce a strong coupling 
to the condensate if the superconductor elastic mean free 
path is sufflciently small. However, disorder in our sim- 
ulations is most likely not strong enough to correspond 
to this possibility. Finally, the Thouless energy of a nor- 
mal cavity appears also in circuit theory^-. Our model 



for the conductance fluctuations clarifies the concept of 
Thouless energy intrinsic to a superconductor. 
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